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Q . Abstract 

■ A nonlocal covariant extension of the two-flavor Nambu and Jona-Lasinio (NJL) model is 
constructed, with built-in constraints from the running coupling of QCD at high-momentum 
and instanton physics at low-momentum scales. Chiral low-energy theorems and basic cur- 
rent algebra relations involving pion properties are shown to be reproduced. The momentum- 
dependent dynamical quark mass derived from this approach is in agreement with results 

Oh| from Dyson- Schwinger equations and lattice QCD. At finite temperature, inclusion of the 

Polyakov loop and its gauge invariant coupling to quarks reproduces the dynamical entangle- 
ment of the chiral and dcconfinement crossover transitions as in the (local) PNJL model, but 
now without the requirement of introducing an artificial momentum cutoff. Steps beyond the 
I mean-field approximation are made including mesonic correlations through quark-antiquark 

■ ring summations. Various quantities of interest (pressure, energy density, speed of sound 
I etc.) are calculated and discussed in comparison with lattice QCD thermodynamics at zero 

■ chemical potential. The extension to finite quark chemical potential and the phase diagram 
in the (T, /i)-plane are also discussed. 

00 ! 1 Introduction 

O 

^ . The thermodynamics of strongly interacting quark-gluon matter is a topic of fundamental in- 



X 



terest from several perspectives, including the physics of the early universe and ultra-relativistic 
^ ' heavy-ion collisions. Modeling the phases of QCD in comparison with results from lattice sim- 

^ ■ ulations has become an active part of this field. A promising approach in this direction is the 

synthesis of the time-honored Nambu and Jona-Lasinio (NJL) model [1,2] with the Polyakov 
loop [3]. This so-called PNJL model [4-11] is based on two principal symmetries: the chi- 
ral SU(A'^f)L X SU(A'^f)R symmetry of QCD with A'f massless quark flavors, and the symmetry 
associated with the center Z{3) of the color gauge group. 

The NJL model offers a schematic but nonetheless quite realistic picture of the basic dynamics 
behind spontaneous chiral symmetry breaking, with the chiral (quark) condensate and pions as 
Goldstone bosons emerging from a chiral invariant, local four-point interaction between quarks. 
Confinement is absent in the NJL model. The Polyakov loop, on the other hand, is introduced as 
the order parameter describing the confinement-deconfinement transition in pure-gauge QCD. 
In the PNJL model, the minimal gauge invariant coupling of the Polyakov loop to quarks 
produces a characteristic dynamical entanglement of the chiral and dcconfinement transitions. 
This phenomenon is indeed observed in recent lattice QCD computations [12], while at the same 
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time its more detailed assessment is still under dispute [13]. Irrespective of such remaining open 
questions, the quahtative comparison of a variety of PNJL model results with those from lattice 
QCD has so far turned out to be surprisingly successful. 

The price paid for the simplicity of the NJL model, with its local point couplings between 
quarks, is the necessity of introducing a momentum cutoff to regularize loop integrals. The 
model draws a schematic picture of the QCD interaction by reducing it to a constant at low 
momenta, \p\ < Anjl ~ 0.6-0.7 GeV, and simply turning off the interaction for \p\ > Anjl- The 
applicability of the model is therefore restricted to energy and momentum scales (temperatures, 
chemical potentials) small compared to Anjl- Connections with the running QCD coupling and 
the established high-momentum, high-temperature behavior governed by perturbative QCD are 
ruled out right from the start. 

In the present work we develop a nonlocal generalization of the PNJL model with running 
coupling, capable of describing both chiral and deconfinement transitions, but without the cutoff 
of the standard NJL model. This is accomplished by making contact with the Dyson-Schwinger 
approach to QCD [14]. Nonlocal extensions of the NJL model are designed to remove deficiencies 
of the local theory, while at the same time the nonlocal interactions regularize the model in such 
a way that anomalies are preserved [15], charges are properly quantized and currents conserved. 
The effective interaction is finite to all orders in the loop expansion and, hence, there is no 
need to introduce extra cutoffs. Our approach takes over those benefits; it differs from other 
works (see e. g. Refs. [16-20]) in that we do not use a separable interaction but derive an explicit 
form for the effective interaction guided by full QCD. The running coupling is constructed in 
accordance with perturbative QCD at high-momentum scales. The interpolation down to low 
momenta is performed in close correspondence with the instanton model. Finally, the Polyakov 
loop is implemented as in the local PNJL model, through the introduction of an additional 
homogeneous temporal background field with standard SU(3)c gauge invariant coupling to the 
quarks. 

This paper is organized as follows. In Section 2 we first construct a nonlocal extension of the 
NJL model at zero temperature. A detailed derivation is given of a nonlocal effective interaction 
constrained by Dyson-Schwinger and lattice QCD results. The momentum dependent dynamical 
quark mass is calculated. Pion properties are derived, and it is demonstrated that fundamental 
chiral symmetry and current algebra relations are properly reproduced. The thermodynamics 
(restricted to vanishing chemical potentials) and the implementation of the Polyakov loop is 
described in Section 3. The resulting nonlocal PNJL model is applied to calculate the chiral and 
deconfinement crossover transitions. The outcome is discussed in comparison with lattice QCD 
results and previous calculations. Pressure and energy density, conformal measure and speed 
of sound are computed around the transition temperature, first in mean-field approximation 
and then with inclusion of mesonic modes (primarily the pion). The extension to finite quark 
chemical potential is briefly outlined. In Section 4 we present our conclusions and an outlook. 

2 Nonlocal Nambu— Jona-Lasinio Model 
2.1 Nonlocal NJL action constrained by QCD 

Consider as a starting point the basic quark color current, Ja = V'(x)7^^V(a;), with the quark 
fields il^{x) and the generators Aa/2 of the SU(3) color gauge group. The QCD effective inter- 
action between two such color currents contributes a term of the following generic form to the 
action: 



where G'^ti^ ~ v) is proportional to the gluonic field correlator. Approximate solutions for this 
correlator can be obtained from Dyson-Schwinger calculations [14] and from lattice QCD [21]. 
We do not need to specify Q^y in all detail. We just note that Eq. (2.1) defines a nonlocal effective 




(2.1) 
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interaction between quarks. The range of the nonlocality is determined by the correlation length 
characteristic of the color exchange through gluon fields. 

Consider next the diagonal parts of Q^^, that have a leading behavior Q1^y{x—y) ~ g^y 6abG{x— 
y), with a distribution function Q{x — y) of mass dimension 2. Then perform a Fierz exchange 
transformation which turns iSint into a sum of terms 

Sint = ^ca [ d^x I d''y,p{x)ra^P{y)g{x-y)i;{y)r^^{x), (2.2) 

with well-defined Fierz expansion coefficients Cq,. The are a set of Dirac, flavor and color 
matrices, resulting from the Fierz transform, with the property 7o-fa7o = r^- By a change of 
variables, Eq. (2.2) can be rewritten as: 

<Sint = ^C„ j d^X j dSV;(x+|)r,^(x-0 ^(z)V;(x-0r«^(x + |). (2.2') 

Throughout this work we restrict ourselves to the two-flavor case with i\){x) = {u{x) , d{x))~^ 
and the quark mass matrix rh = diag(m„, m^^). We will work in the isospin-symmetric limit 
with ruu = rud = ruq. The Fierz transformation introduces combinations of operators -Tq which 
maintain the symmetries of the original QCD Lagrangian. The symmetry that governs low- 
energy QCD with two flavors is global chiral SU(2)l x SU(2)r,. A minimal subset of operators 
satisfying this symmetry is the color-singlet pair of scalar-isoscalar and pseudoscalar-isovector 
operators, = (1,175-?), the one we focus on. Other less relevant operators (vector and 
axialvector terms in color singlet and color octet channels) will be ignored in the present paper. 

Our nonlocal NJL model is thus defined by the following Lagrangian: 

C{x) = i^{x){i0 - m)iP{x) + Cintix) , (2.3) 
Ant(a;) = J d'^zg{z) [VS(x+)V'(a;_)V5(x_)V'(a;+) -|-'i/;(x+)i75fV'(x_) • 'ip{x-) 1-^5 Tip {x+)] , 

with x± = X ± z/2. The action in Minkowski space with x^ = {xo,x)~^ is S = J d^xC From 
here on and throughout this work calculations will be performed in Euclidean space with the 
corresponding action 

5e = y" d^xC^ = j d^x^{x){-i$ + m)^{x)+St^ , (2.4) 

where x = (.■?, X4) etc. are now Euclidean space-time coordinates. For convenience we replace 
the function Q{z) representing the nonlocal interaction kernel by a coupling constant times a 
normalized distribution, 

g{z) = ^C{z), (2.4a) 

with J d'^zC{z) = 1. The coupling strength G of dimension (length)^ (i.e. mass dimension —2) 
is understood to absorb the corresponding Fierz coefficient. Note that the standard (local) NJL 
model follows for the limiting case C{z) = 5'^{z). 

The bosonization of the Euclidean action is then performed as usual by introducing aux- 
iliary fields (pa{x) = {a{x),TT{x)), where a and vf are chiral partner boson fields representing 
scalar-isoscalar and pseudoscalar-isovector mesonic degrees of freedom, respectively. Writing 
the generating functional in path integral formalism, Z = J ^-i/^^ip exp[— (Se], we finally obtain 

Z = J ^a^Tf exp[-5^°^] , (2.5) 

with 

= - In det 1 + ^ 1 ^ ip) ip) , (2.6) 
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where we have introduced 

Eq. (2.6) involves the momentum space representation of the operator A, as fohows: 

Aip,p') := {p\A\p') = i-p + mg) (27r)4,5(4)(p-p') T^Re [4>a{p-p')] , (2-7) 

with 



C{p) = J dS e-'P-'^Ciz 



We use the notation C{p) =: C{p) for simphcity and note that C{p = 0) = 1. Details of the 
bosonization procedure are given in Appendix A. 

2.2 Mean field approximation and beyond 

Next, assume that in the homogeneous and isotropic vacuum, the scalar a field has a nonzero 
expectation value a = {a), while the vacuum expectation values of the pseudoscalar fields ttj are 
zero. We write cr{x) = a + 5a{x), Tf{x) = 6tt{x) and expand the bosonized action (2.6) around 
the mean field in powers of the mesonic fluctuations 6a, Sir: 

S^°'=S^^ + S^^ + ... (2.8) 
The mean field contribution per four-volume V^'^^ is given by 

cMF r j4 -2 

with the mass function M{p) determined by the gap equation 

M{p) = mq + C{p) a . (2.9a) 

The quadratic terms beyond mean field approximation are derived explicitly in Appendix B. 
Here we only state the result: 

where 

.±/.2^_ 1 /■ d'*? , . q-{q + p)^M{q)M{q+p) 



F^{p') = ^-8N,J A^C{q)C{q+p) 



[q2 + M^q)][{q+py + M^q+p)] 
Q 




q + p 



(2.11) 



and r+jP denote the scalar-isoscalar and pseudoscalar-isovector vertices, respectively. The 
loop diagram involves the fermion (quark) quasiparticle propagators (in Euclidean space^) 

Sf(p)= = Kttt (2-12) 

with the momentum-dependent (dynamical) constituent quark mass M{p) (cf. Sect. 2.6). 



^We use the convention p4 := ipo, 74 := i7o- 
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2.3 Gap equation and chiral condensate 

The mean field part, S^^, of the action Se is governed by the scalar mean field a. Its value is 



found by the principle of least action, ^ = for cr = a. It is straightforward to derive 



a = 8NcG 



, M{p) 



(2.13) 



to be solved self-consistently with the gap equation (2.9a). The chiral condensate (V'V') = 
{uu) + {dd) can be calculated from 5^°^ using the Feynman-Hellmann theorem, by differentiation 
with respect to the current quark mass ruq. 
Equivalently, one can use the definition 



-iTr lim 



with the full fermion Green function. 



(2.14) 



(2.15) 



and the free quark propagator Sp^ subtracted. Using Eq. (2.12) this leads to 



J (27r)4 [p2 + M2(p) p'^ + ml 



M{p) 



(2.16) 



Note that M{p) rUq for large p. The subtraction makes sure that no perturbative artifacts 
are left over in (V-'V') foi" '^g 7^ 0. 

At this point a comparison with the local NJL model is instructive. There, C{p) and M(j)) 
are replaced by constants, C(0) = 1 and M. The relevant momentum space integrations are cut 
off at \p \ = Anjl- The gap equation is simply M = —G{'tpilj) = a. This direct proportionality 
between the scalar mean field a and the chiral condensate is not realized any more in the nonlocal 
model, as the inspection of Eqs. (2.13) and (2.16) demonstrates. 

2.4 Pion mass, queirk-pion coupling constant and pion decay constant 

The pion mass is determined by the pole of the pion propagator, while the square of the 
quark-pion coupling constant, g^qq, figures as the residue of the pion pole in the pseudoscalar- 
isovector quark-antiquark amplitude, 



15 



Trqq 



+ m. 



2 175Ti . 



(2.17) 



Its inverse is easily calculated from the bosonized action (2.6) by means of functional differen- 
tiation, i. e. 

The last equality results from Eqs. (2.10) and (2.11). The pion mass is therefore determined by 

F-{-ml)=Q. (2.19) 
Furthermore, by comparison with Eq. (2.17) the quark-pion coupling constant is given as: 

dF-{q^) 



-2 
dnqq 



dg2 



(2.20) 



5 



The calculation of the pion decay constant /tt is more involved. It is defined through the 



matrix element of the axial current J^{x) between the vacuum and the physical one-pion state 

TTj with T^j = g-^qT^j-. 



(0|J^(0)|7f,(p))=i5,,p^A. (2.21) 

Note that the usual local axial current, Jj^''°'^ = 'i/'(rc)75^7^^/'(a:), is not conserved within the 
nonlocal framework (cf. Ref. [22]). In order to obtain an explicit expression for the axial current, 
we gauge the action <Se by introducing a set of axial gauge fields A^^{x) in the standard way; 
partial derivatives arc replaced by covariant derivatives, i.e. ^ + ^75 Tj^^x), and due to 
the nonlocal nature of the theory, interaction currents have to be replaced by [17, 23, 24] 



with the function 



W{x,y) = V jexp 



(2.22) 



1 
2 

including the path-ordering operator V. This assures that the gauged action Sj^ is invariant 
with respect to local axial transformations. 

After those replacements a gauged bosonized action, 5^°**, can be built in the same fashion 
as outlined in Appendix A. Prom there it is easy to calculate the matrix element required for 
the determination of /^r according to Eq. (2.21): 



r2 cbos 



A=0 



Finally, switching to momentum space and evaluating the functional derivatives and traces 
according to Appendix B leads to the pion decay constant in mean field approximation^ (see also 
Ref. [17]): 

/tt = 9Tvqq 2 ' (2.23) 



where we have defined 



with 



J (27r 



Xi = 8Ar, / —^C{q) 



'^Q ni \ q- (q + p) + M{q)M{q+p) 



)4 ^■^^[g2 + M2(g)] [(g + p)2 + M2(g + p)] 
d^g , X q- (q+p) + M{q)M{q+p) 



l2 = m,j ±^^C{q)C{q+p)- 



{2ttY ^^[g2 + M2(g)][(g+p)2 + M2(g+p)] • 

Two pion observables in vacuum, and mjr, together with g^^gg will be used to fix the parameters 
of the nonlocal model of this work (see Sect. 2.7). 



2.5 Chiral relations 

In order to confirm that our model is consistent with fundamental chiral symmetry requirements, 
we derive the Goldberger-Treiman and Gell-Mann-Oakes-Renner relations. For this purpose 
it is sufficient to consider an expansion of the integrals Ii,l2 up to 0{m1), i.e. to leading 
(linear) order in the quark mass rriq. First note that using Eqs. (2.13) and (2.19), one can write 
jr(0) = a/G and l2{-ml) = 1/G. Then one obtains from Eq. (2.23) 

mlU = rUg g^^qq Ii {-ml). (2.24) 

^Note, that the result stated here originates from an RPA calculation with vertex functions chosen such that 
Ward identities are preserved (see Ref. [22]). The complete result, however, is more involved and its derivation is 
relegated to Appendix C. The result (2.23) corresponds to the leading order expansion in rrin of (C.3). 
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A Taylor expansion about the chiral limit of the pion polarization term, Eq. (2.11), leads to 

dF- 



F-{p') = F-{0) + 



+ Oip"") = F- (0) + g-g\ p' + 0(/) 

p2=0 



Using this in the conditional equation for the pion mass, F (— m^) = 0, one recovers the leading 
order contribution to mi: 



"IT- 
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mi = ^m,Xi + 0{mi) = F'iO) + 0{mt) . (2.25) 



a 



Finally, combining Eqs. (2.11) and (2.16) gives 

F-iO) = -'^m+Oiml). 
Prom Eq. (2.24) and the first equality of Eq. (2.25) one finds 

Ugnqq = ^ + Oimg), (2.26) 

which, in the chiral limit, i. e. for mq = 0, is nothing but the Goldbergcr-Trciman relation at the 
level of quarks as quasiparticles. Using this relation, the second equality of Eq. (2.25) together 
with the expression for F^ (0) derived above gives 

mlf^ = -mq{i,i^)+Oiml), (2.27) 

which is the well-known Gell-Mann-Oakes-Renner relation. We have thus demonstrated that 
the nonlocal NJL model preserves chiral low-energy theorems and current algebra relations. 



2.6 Determination and fixing of the distribution C{p) 

The Fourier transform C{p) of the nonlocality distribution introduced in section 2.1 is a key 
quantity and the basic input of the present approach. QCD constraints are used to determine 
C{p), both in the high-p (perturbative) and low-p (non-perturbative) regions, as follows. 



2.6.1 Quark self-energy at large momentum 

Consider a quark propagating with large (Euclidean) momentum p in the QCD vacuum. Its 
self-energy I^{p) is given pictorially as 



-iU{p) = 




(2.28) 



with full quark and gluon propagators and vertex functions^. At high momentum, this self- 
energy can be evaluated recalling its operator product expansion and identifying the leading 
0{p-^) term [25,26]. 

It is well known that QCD perturbation theory at any order does not generate a non-zero 
quark mass term starting from a massless quark, simply as a consequence of helicity conservation 
at each quark-gluon vertex. However, as pointed out by Politzer [25], the presence of a non- 
perturbative QCD vacuum with a quark condensate {'ipip) ^ turns a massless current quark 
into a quasiparticle with a momentum dependent mass (often referred to as a constituent quark): 
the spontaneous breaking of chiral symmetry implies a non-zero dynamical quark mass. 

^In a self-consistent Schwinger-Dyson approach the filled circles denote fully dressed propagators and the open 
circle is the irreducible quark-gluon vertex. 
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At high quark momentum, the leading non-trivial part of ^{p) is found replacing the full 
quark-gluon vertex function by the bare vertex ( "rainbow truncation" ) and replacing the full 
gluon propagator by its free part, i.e. 4TTas{p'^)Dfi^{p) — > 4-Kasip'^)D^^^^'^{p) ("Abelian approxi- 
mation"), with the running QCD coupling as(p^)- The full quark propagator in (2.28) is kept 
with its non-perturbative pieces. One then arrives at the quark self-energy to leading order in 
(withp^ > 0): 

= 2^^^^ + ^^^^^^ + ' ^^-^^^ 

where we have explicitly written the non-perturbative piece, S'^' P' , and the remaining 6S stands 
for all perturbative corrections to E (which vanish in the chiral limit, ruq ^ 0). The parameter 

^ is required for gauge fixing (see also Refs. [26, 27]). 

From Eq. (2.29) one arrives at the dynamically generated quark mass term 

M{p') := E^-^-ip') = p2 + ■ (2.30) 

This result coincides with the constituent quark mass calculated by Politzer. As mentioned, even 
in the limit of vanishing current quark masses, the quark propagator has a non-perturbative 
contribution which generates a large constituent quark mass. In what follows we choose the 
Landau gauge, setting ^ = 0. Any other value of ^ can be absorbed by a rescaled coupling in 
Eq. (2.30). 

From Eqs. (2.4a) and (2.12) it is evident that the distribution C{p) must behave as M{p) at 
large p, up to normalization. In order to reproduce the QCD result (2.30) in the chiral limit 
(with Nc = 3 and N{ = 2), we therefore choose 



C(p=)oc^^^ (P>r), (2.31) 

6 p 



at p larger than a matching scale F of order 1 GeV below which C{p'^) is governed by non- 
perturbative physics (see the following subsection). For our present purposes it is sufficient to 
use the standard form^ 

4-7r 

«s(p') = ^ , o (2.32) 

with /3o = 9 (using three active quark flavors) and Aqcd = 0.25 GeV in order to reproduce 
as = 0.12 at mz = 91.2 GeV. 



2.6.2 Aspects of non-perturbative QCD 

The behavior of C(p) at low momentum p, below the matching scale F, is not controlled by first 
principles and needs guidance e. g. from lattice QCD and other non-perturbative considerations. 
A useful scheme for extrapolating to low momenta is provided by the Instanton Model (IM) 
[28]. Instantons are non-perturbative gauge configurations representing extrema of the Yang- 
Mills action. They are known to play an important role in the Euclidean partition function. In 
particular, it has been demonstrated that instantons provide a mechanism for chiral symmetry 
breaking in the QCD vacuum. Detailed computations in the context of the Interacting Instanton 
Liquid Model [29] have shown that, in the chiral regime of QCD, there is not only qualitative but 
also quantitative agreement between the IM and full QCD, at least for = 2. It is also known 
[28] that the partition function of the IM can be represented in terms of a nonlocal fermionic 
effective action similar to that of the nonlocal NJL model, where now the nonlocality originates 
in the finite size of the instantons (~ 1 /3 fm) which mediate the interaction. 

The more precise NNLO form for as(p^) can of course be used but the resulting changes are insignificant 
within the accuracy of the present approach. 
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The distribution Ci{p) extracted from the IM is expressed in terms of Bessel functions as 

[28] 

Clip) = 7rp''d^^iIo{x)Koix) - hix)Ki{x)) , ^ = ^ , (2.33) 

where d ~ 0.35 fm is the typical instanton size. The implementation of this instanton inspired 
Cj{p) in finitc-tcmperature calculations is numerically expensive. For practical convenience we 
approximate C(j)) at low momentum p by a Gaussian, 



C{p^) = exp 



p^d 



p<r, (2.34) 



subject to the requirement that this C{p) smoothly continues into the logarithmic high- momentum 
behavior at the matching scale F. The width of the Gaussian {d~^ ~ 0.56 GeV) is then directly 
determined by the instanton size. 

Figure 1 shows a comparison of the Gaussian form (2.34) and the IM result (2.33) for C{p), 
demonstrating that the two curves are indeed quite similar. One should note in addition that, 
when implementing C{p) in the gap equation (2.13) determining the dynamical quark mass 
M{p), or in the calculation of the chiral condensate (2.16), the results are insensitive to the 
detailed behavior of C{p) or M{p) at very low momentum. The reason is that the Euclidean 
momentum space integrals (2.13) and (2.16) introduce a measure proportional to dpp^ so that 
the maximum weight under those integrals is concentrated in a characteristic window around 
p - 0.8 GeV. 

2.7 Pcirameter fixing and dynamical qucirk mass M{p) 

Given the nonlocality distribution 

U-p'dy2 forp2<r2 

C{P')={ ^ a^ip') , 2 (2.35) 
const. • — for P > T'' , 
K p^ 

normalized as C(j) = 0) = 1 and with a constant fixed by the matching condition at p = F, 
there remain only two parameters apart from the matching scale F: the coupling strength G in 
Eq. (2.4a) and the (current) quark mass mq as a measure of explicit chiral symmetry breaking. 

The prime advantage of the present nonlocal approach is the absence of an artificial momen- 
tum space cutoff as it appears in local NJL type models. Indeed, the gap equation (2.13) does 
not require any regularization: the scalar mean field a and, likewise, the dynamical quark mass 
M{p) are well-defined within the nonlocal scheme.^ 

The parameters and resulting values of observables are listed in Table 1. The determination 
of G and ruq is such as to reproduce, as closely as possible, the empirical pion decay constant 
/tt = 92.4 MeV and the pion mass mT^. Fig. 2 shows the momentum dependence of the dynamical 
quark mass M{p) compared with lattice data from Ref. [30]. One notes that the nonlocal NJL 
model reproduces very precisely both the low- and the high-momentum behavior. 

The optimal matching scale F ~ 0.8 GeV turns out to be located in the momentum window 
where the gap equation (2.13) has its maximum weight. It is therefore important to examine 
the sensitivity of the results with respect to variations of F. We have performed this test by 
varying the matching scale in the range 0.6 GeV < F < 1 GeV (i. e. about 20 % around its 



■''This docs not exclude the possibiUty that secondary quantities, such as the chiral condensate, can be weakly 
divergent and potentially require a cutoff at ultrahigh momenta, very fax beyond the range of applicability of the 
model. Inspection of the integral for (ijip) in Eq. (2.16), with the asymptotic form M{p) — > mq + ^i^^^^i inserted, 
displays a weak (double-logarithmic) far-ultraviolct divergence. For regularization we choose a cutoff at 20 GeV, 
the necessity of which just reflects the simple leading-order choice for the asymptotics of C{p). 
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— non-local NJL 

- Instanton model 
■ ■ local NJL 



1.0 1.5 2.0 

p [GeV] 



2.5 




1.0 1.5 
p [GeV] 



Figure 1: Comparison of the form factor C 
used in the nonlocal NJL model (solid line) 
and the analogous form factor obtained from 
the Instanton Model (dashed line). The dot- 
dashed line shows the function C{p) in an NJL 
model with four-momentum cutoff. 



Figure 2: Momentum dependence of the con- 
stituent quark mass M compared with lattice 
data extrapolated to the chiral limit (from 
Rcf. [30]). The dashed line shows the mass 
M calculated fully within the IM for the chi- 
ral hmit [28]. 



ruq [MeV] 


r [GeV] 


G [GeV-^] 


3.3 


0.83 


41.1 



(V'V)^^^ [GeV] 


a [GeV] 


U [GeV] 


[GeV] 


-0.36 


0.42 


0.09 


0.14 



Table 1: Parameters and calculated physical quantities for the nonlocal NJL model with Nc = S 
and iVf = 2. 



optimal value), with the constraint of keeping the chiral condensate (V'V') fixed. The resulting 
variations in the pion mass and decay constant are within only 10%, implying stable conditions. 

Note that a direct comparison to the NJL model is not appropriate, because the NJL model 
uses a three-momentum cutoff for the regularization of the integrals while in the nonlocal ap- 
proach the distribution C{p) mimics a four-momentum cutoff. This implies, in particular, a 
different scale for the coupling strength G appearing in the scalar mean field a. Consider 
again the chiral limit and compare Eqs. (2.13) and (2.16). One can write the dynamical 
quark mass as M{p) = C{p)a = -C{p)CG{i)i}) where C = - J d'^pC{p)f{p)/ J d'^p f{p) with 
f{p) = M{p)/ (p"^ + M'^{p)). Given that the maximum weight under the integrals is located 
around p ~ 0.8 GeV, one finds C ^ ^ and, consequently, M(0) ~ —^{iptp). Hence with the 
nonlocality profile as shown in Fig. 1, the coupling G in the present approach should be about 
four times that of the standard NJL model, which is indeed the case. 

The quark mass 2± 3.3 MeV is compatible with QCD estimates at a typical renormal- 
ization scale n ~ 2 GeV [31]. The value of the chiral condensate follows consistently from the 
Gell-Mann, Oakes, Renner relation (2.27), with {uu) = {dd) ~ -(0.28 GcV)^. The matching 
scale r ~ 0.85 GeV for the low- and high-momentum parts of C(p) can be considered "natural". 



3 Thermodynamics of the nonlocal NJL Model 

In this section we apply the nonlocal NJL approach in order to model QCD thermodynamics. 
Of primary interest are the chiral and confinement-deconfinement transitions. This requires cou- 
pling the quark quasiparticles to the color gauge fields. One-gluon exchange degrees of freedom 
(at the level of the quark self-energy) are effectively incorporated in the high-momentum part 
of C{p), while chiral dynamics is implemented through the scalar field a which plays the role of 
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a chiral order parameter. The thermodynamics of deconfinement, on the other hand, is intro- 
duced by couphng the quarks to the Polyakov loop as the order parameter of the confinement- 
deconfinement phase transition in pure-gauge QCD. This is the step to the Polyakov-loop- 
extended NJL (or PNJL) model [4-11]. We now proceed to develop the nonlocal generalization 
of this PNJL model. 

3.1 Polyeikov loop 

The so-called renormalized Polyakov loop, ($), is defined through 

(^.(f)) = i-(tre[L(x)]), (3.1) 

where trc denotes the trace over color space only and L is the Polyakov loop, a timelike Wilson 
line connecting two points at i = and t = —i(3 in imaginary time space, with periodic boundary 
conditions: ^ 

L{x) =Vexp!^i J dryl4(r,f)| . (3.2) 

Here V is the path-ordering operator and A4 := A'^ta = i^o^a is the fourth component (in 
Euclidean space) of the gluon fields^. 

As shown e.g. in Ref. [32], ($) = implies an infinite free energy, corresponding to confine- 
ment, while ($) = 1 implies a vanishing free energy of the system and, thus, deconfinement. Ac- 
cording to Ref. [33], ($) serves strictly as an order parameter for the confinement-deconfinement 
phase transition only in the case where no (not even massless) quarks are present. In other cases 
($) is still useful as an indicator for a rapid crossover transition. 



3.2 Coupling of queirks and PolyEikov loop 

Gluons can be introduced into the nonlocal NJL model by the standard minimal gauge coupling 
procedure, substituting the normal derivate by the covariant derivative^ = — iA^ or, 
in momentum space, + A^.^ 

As usual we introduce explicitly particle and antiparticle degrees of freedom switching to 
Nambu-Gor'kov space. This is particularly useful for the correct gauging of the theory. The 
fermion propagator is derived from the Euclidean action 5e, Eq. (2.4), using standard functional 
calculus. The inverse propagator is given by the fermion determinant already calculated in 
Eq. (2.7). Hence, in mean field approximation we have 

SmM = -t + mq + C{jp)a . (3.3) 

The transition to Nambu-Gor'kov space consists in the introduction of the field := , 
where tp'-^ = C'lp'^ with the charge conjugation matrix C = 17072, and the identification 




with r = ^c),r^ = -crc. 

"The SU(3) matrices ta,a € {1, ...,8} are defined through the Gell-Mann matrices Aa,o € {1,...,8} as 

t — ^ 

'^The coupUng strength g is assumed to be included aheady in the definition of the gluon fields, i. e. gA^ — » A^. 

*This minimal substitution procedure can be derived rigorously as follows: the (gluon) gauge fields are intro- 
duced in the nonlocal model according to Eq. (2.22), with 75r''.4Jl replaced by in the exponential. Limiting 
ourselves to the case of constant gauge fields A4, Af, the path- integrations in the exponential can be carried out, 
generating terms e'^'''**. In momentum space, these terms translate into delta functions 6{qA — {pA — A\)) S^^\q—p) 
leading to the substitution prescription appropriate for the nonlocal approach. 
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The introduction of the Nambu-Gor'kov propagator is particularly convenient when finite 
temperatures are considered. After bosonization of the Lagrangian the partition function of the 
theory becomes (with (3 = 1/T the inverse temperature): 



In detA 



exp 



(3.5) 



where ipa represents the mesonic degrees of freedom and A = 5*^ g^^^g^ is the Fermi propagator 
coupled to the gauge fields, obtained from Eq. (2.7) after the substitution Pfi + A^j,. 

Without loss of generality the Polyakov loop field can be represented in terms of the a = 3 
and a = 8 color components of the gluon fields, i. e. we consider only the diagonal matrices of the 
SUc(3) Lie algebra (the so-called Cartan sub-algebra of SUc(3)) and neglect spatial fluctuations 
at this stage so that only the time component is relevant: 

A^ = 5^i{Alh + A%). 
The group integration over the remaining non-diagonal elements gives the Haar volume 



1 



SU(3) 



n = (cOs(^^3) - COs(\/308)) 

ie{l,2,4,5,6,7} 



sm 



(3.6) 



where we have set (f)s,8 = This volume can be written in terms of the Polyakov loop^ $ 

and its complex conjugate $*, 
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with $ = -jy-trc [exp (i(^3A3 + ^s-^s))]- With the partition function (3.5) this procedure leads 
to the construction of the following effective potential U that includes the effects of the six 
non-diagonal gluon fields: 



J($, $*) = _^ 1 _ 6$*$ + 4 (|$*^ + - 3 ($*$)^ 



(3.7) 



m:^^ = -l62(r)$*$ + 64(T) In 
The coefficients are parametrized as 



1 - 6$*$ + 4 (^^*^ + - 3{^*<^f 



(3.8) 



62(r) = ao + ai ( ^ j + 02 f ^ 



hiT) 



+ as 



The first term on the right-hand side is reminiscent of a Ginzburg-Landau approach. The values 
of the coefficients are taken from Ref. [8] and listed in Table 2. The parametrization (3.8) of the 
Polyakov loop effective potential U is applicable at temperatures T up to about twice the critical 
Tc- At higher temperatures, transverse gluon degrees of freedom - not covered by the Polyakov 
loop - begin to be important. Other parametrizations of U are possible. An example is the 
two-parameter ansatz [34] based on the strong coupling limit. These two versions of U differ at 
high temperatures but produce very similar pressure profiles [34] at T < 2Tc, the temperature 
region of primary interest in the present study. 

We now continue with the evaluation of the partition function (3.5). In the following we first 
limit ourselves to the mean field approximation. The integrals in Eq. (3.5) are easily performed. 
The thermodynamic potential J7 = — ^ In Z becomes 



n 



■ In det 



(3.9) 



'EYom now on we omit angled brackets for notational simplicity, i. e. we write ($) —^ 3>. 
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ai 


02 


03 


64 


3.51 


-2.5G 


15.2 


-0.62 


-1.68 



Table 2: Parameters of the Polyakov potential U (from Ref. [8]). 



where S^^ is the inverse quark propagator expressed in Nambu-Gor'kov space and given in 
Eq. (3.12). The functional determinant is evaluated using the Matsubara imaginary time 
formalism (sec e.g. Ref. [35]), replacing pQ — > icun, with the fermionic Matsubara frequencies 
LOn = (2n + 1)ttT, (n G !Z), and performing integrals over four-momentum space according to 



(zero temperature) 



d^j9 

(27r)4 



(2^ 



(finite temperature) . (3.10) 



The resulting thermodynamic potential is^°: 

d^p 



(2t)- 



■ tr In 



(3S-\iujn,p)\+^+U{<^,<P*,T), 



(3.11) 



with 



S ^{iujn,p) = 



it^nTo - J ■ p- M - 1^4470 







ia;„7o - 7 M* + 1^4470 

where the momentum dependent mass matrix M is diagonal in color space, 

M = dmg,{M{uj-,p),M{u+,p),M{iolp)) , 
with = LOn i A4,uj^ = <^n- The trace may be further simplified leading to 



i=0,±neZ 



d^p 
(2^ 



In 



col, +r + M\uj\,p) +_+Z^($,$*,r) 



(3.12) 



(3.11') 



This is the thermodynamic potential of the nonlocal PNJL model in mean field approximation^^. 



3.3 Gap equations and results 

Apart from the sigma field, the nonlocal PNJL model includes two more degrees of freedom, 
$ and or ^ and A|, that mimic the thermodynamics of confinement and deconfinement. 
Therefore, three gap equations are derived and solved for a, # and The necessary conditions 
are given by the requirement of a stationary potential: 



89. dn 



on 



da dAl dAl 



0. 



(3.13) 



Following Refs. [8,10], we have $ = $* in the mean field approximation and, consequently, 
^1 = 0. 

In Fig. 3 we show the results for the temperature dependence of the chiral condensate and of 
the Polyakov loop using the parameters given in Tables 1, 2. This figure illustrates once more, 

as already demonstrated in Refs. [8,10], the entanglement of chiral dynamics and Polyakov 
loop degrees of freedom, a characteristic feature of the PNJL approach. In the absence of a 



^''Note an extra factor i because of the doubling of the degrees of freedom in Nambu-Gor'kov space. 
^^Diquark degrees of freedom are omitted in the present work and will be dealt with elsewhere. 
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Figure 3: Solid curves: calculated temperature dependence of the chiral condensate (V'V') and 
of the Polyakov loop $ normalized to the critical temperature Tc = 207 MeV as obtained in the 
nonlocal PNJL model. Also shown for orientation arc lattice results from Ref. [12]. The dashed 
lines show the chiral condensate in the chiral limit (left for the pure fermionic case) and of the 
Polyakov loop for the pure gluonic case, respectively. 



coupling between quark quasiparticlcs and Polyakov loop, and in the chiral limit, the second 
order chiral phase transition (for A^f = 2 flavors) and the first order deconfinement transition (of 
pure gauge QCD) appear at very different critical temperatures (Tchirai ~ 110 MeV for the chiral 
phase transition and Tq 270 MeV for deconfinement). The presence of quarks breaks the Z(3) 
symmetry explicitly and turns the first order deconfinement phase transition into a continuous 
crossover. The quark coupling to the Polyakov loop moves this transition to lower temperature. 
At the same time the chiral transition (with explicit symmetry breaking by non-zero quark mass) 
turns into a crossover at an upward-shifted temperature, just so that both transitions nearly 
coincide at a common temperature Tc ~ 200 MeV. 

This symmetry breaking pattern seems also to be realized in recent lattice QCD results 
[12] where a common chiral and deconfinement transition temperature = (196 zt 3) MeV is 
observed. Fig. 3 shows these lattice data for orientation^^. 

3.4 Pressure and mesonic corrections 

So far the calculations have been performed in the mean field approximation in which the 
pressure P = — f7 is determined by the quarks moving as quasiparticlcs in the background 
provided by the expectation values of the sigma field, a, and the Polyakov loop In order to 
get a realistic description of the hadronic phase (at temperatures T < Tg), it is important to 
include mesonic correlations. The hadronic phase is dominated by pions as Goldstone bosons 
and their interactions. The quark-antiquark continuum is suppressed by confinement. However, 
as pointed out in Ref. [36], mesons described within the standard (local) PNJL model can still 
undergo unphysical decays into the quark-antiquark continuum even below Tc- In the nonlocal 
PNJL model, those unphysical decays do not appear by virtue of the momentum dependent 
dynamical quark mass. Following Rcfs. [22] one can easily show that, within our nonlocal 
NJL model, no real poles appear in the fermion propagator (2.12) and hence mesons are stable 

^^A direct compaxison with present nonlocal PNJL calculation is not possible since the lattice computation has 
been performed with iVf = 2 + 1 flavors and a pion mass m-w ^ 200 MeV. 
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well below Tc- This means that the pressure below Tc is basically generated by the pion pole 
with its almost temperature independent position. Therefore the calculated pressure below Tc 
corresponds to that of a boson gas with constant mass. 

To include the mesonic contributions to the pressure in our nonlocal PNJL model we can ba- 
sically use the formalism described in Ref. [36]. One has to calculate the quark loop contribution 
to the mesonic self-energies 777r,o-(i^m)P ) (where = 2TTmT,m G Z is the bosonic Matsubara 
frequency and p is the momentum of the incoming pion or sigma), depicted in Eq. (2.11), at 
finite temperature. Using the rules (3.10) we end up with 

^ coijuji + ^/n^) + k{k +p)± Mjioj + vm,k+ p)M{uji, k ) 

X ~z 



{ui + i^m)^ + {k + pY + M2(< + u^,k+ p)\ [(<)2 + A; 2 + M^{ui, k : 

(3.14) 

where uj^ = a;„ =b ^44, o;*^ = f^n and M{ujn,p) = mq +C{ijJn,p)cr. The additional contribution of 
mesonic quark-antiquark modes to the pressure is given by a ring sum of RPA chains, investigated 
in Ref. [37] and leading to the expression 

Pmeson(T) = -T ^ / 7^ In [1 - G7TM(i^m, #)] , (3.15) 

M=7r,(7 mGZ"^ ^ 

where cLm is the mesonic degeneracy factor, i. e. (It^ = 3, dg- = 1- Due to the momentum depen- 
dence of the nonlocality distribution C (p) and the dynamical quark mass M(p), integrations and 
summations in Eqs. (3.14) and (3.15) can only be carried out numerically. 

Results for the pressure in the presence of pion and sigma mesonic modes are presented 
in Fig. 4. Apart from the full result (solid line) we additionally show the mean-field result 
(with the pressure determined by quark quasiparticles only) and the mean-field result plus pion 
contributions. It is evident that at low temperatures the mean-field contribution stemming 
from the quarks is suppressed and the pressure can be described by a free pion gas. Near 
the critical temperature the sigma mesonic mode gives a small additional contribution. Finally, 
above temperatures T > 1.5 Tc the mesonic contributions become negligible and the quark-gluon 
mean-fields dominate the pressure. 

Finally, in Fig. 5 we show a comparison of the pressure calculated with the physical pion 
mass, m-TT = 140 MeV, and with a "heavy" pion (niT^ = 500 MeV) corresponding to a quark mass 
of order niq ~ 100 MeV that has frequently been used in earlier lattice QCD computations. In 
this case the mesonic contributions to the pressure are evidently reduced. This explains the 
apparent agreement of lattice data with mean field calculations [7, 8] . 

3.5 Related thermodynamic quantities 

Given the partition function (or the thermodynamical potential Q) we can further investigate 
the energy density e, the trace anomaly (e — 3P)/T'^ and the sound velocity Vs (compare also 
Refs. [38]). In particular, the trace anomaly is of interest here since it is the quantity which 
can be directly computed in lattice simulations (Ref. [12]). The trace anomaly, corresponding 
to the trace of the energy-momentum tensor, is the "interaction measure" which - using ther- 
modynamical relations - can be expressed in terms of a derivative of the pressure with respect 
to temperature: 

S-3P j.±(P\ 



dT V / 
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Figure 4: Pressure calculated in the nonlocal PNJL model in units of as a function of 
temperature normalized to the critical temperature. Solid curve: full calculation (i. e. mean 
field result plus mesonic corrections). Dotted curve: mean-field result (no mesonic corrections). 
Dashed curve: mean-field plus pionic modes (no sigma). 



A further interesting quantity is the ratio of pressure and energy density. Finally, the square of 
the sound velocity (at constant entropy S) is deduced as: 



dP 




dP\ 

dT\v 




s 


rp d^P 
J- QfJ 


V 



which is connected to the second derivative of the thermodynamic potential. 

The following figures show the quantities just mentioned as they result in the mean field 
case and with inclusion of mesonic corrections. Again, mesonic corrections are important only 
at temperatures below (Figs. 6-9). Mesonic correlations do have a strong influence on the 
sound velocity below Tc (Fig. 9). At this point it is not yet possible to compare these results to 
lattice data quantitatively since, as pointed out in the previous section, for larger meson masses, 
these corrections become less important and consequently their impact on Vs diminishes. 

3.6 Finite quark chemical potential and phase diagram 

The nonlocal PNJL approach described in this work can be extended to finite quark chemical 
potential We do this here with the aim of drawing a schematic phase diagram in the (T, jj.) 
plane. For the sake of simplicity we restrict ourselves to a scenario without diquarks. Their 
detailed treatment will be reported in a forthcoming paper. The introduction of a chemical 
potential^^ is accomplished using the prescriptions of the Matsubara formalism (see Refs. [35]): 
shift the frequencies Un — *• uJn—ilJ- in the particle sector (i. e. upper-left submatrix) of the Nambu- 
Gor'kov propagator (3.12) and replace a;„ — > a;^ + i/U in the corresponding anti-particle sector 
(i. e. lower-right submatrix). It is then straightforward to compute the thermodynamic potential 

^''For convenience we use a quark chemical potential throughout this work. The corresponding baryon chemical 
potential is three times the quark chemical potential, i. e. hb = 3/Lt. 
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Figure 5: Comparison of the pressure for physical pion mass, = 140 MeV (solid line), and 
heavy pion mass, m-;^ = 500 MeV (dashed line). Note that the crossing of the curves is due to 
the different critical temperatures they were normalized to. 



^{T,fi) at non-zero fi, following Eqs. (3.11), (3.12) with Matsubara frequencies properly shifted 
by the chemical potential. 

We focus here on the T- and /it-dependence of the scalar field a acting as a chiral order 
parameter, deduced from the condition ^^^'^^ = 0. The result is shown in Fig. 10. The profile 
of a displays once again the chiral crossover transition at fi = 0. It turns into a first order phase 
transition at a critical point (here: Tcep = 167 MeV and hcep = 175 MeV). This qualitative 
feature is typical for NJL or PNJL type models with or without explicit diquark degrees of 
freedom (see c. g. Rcfs. [7, 8, 39]). Related work is reported in Refs. [19, 20] where the nonlocality 
of the Fermionic interaction is introduced only in the three-momentum sector. The critical point 
found in these calculations differs both in its T- and /j-value from the aforementioned references. 
Reasons for such discrepancies are given in Refs. [8, 34] where it is pointed out and demonstrated 
that the location of the critical point is extremely sensitive to model details and input parameters. 

Of course, the picture drawn in Fig. 10 is to be taken as only schematic, for several reasons. 
First, the location of the critical point is extremely sensitive to the input current quark mass 
[8]. Moreover, extensions to A^f = 2 -|- 1 flavors including the axial U(1)a anomaly (Rcfs. [34, 40]) 
point out uncercainties not only in the location of the critical point but also for the mere existence 
of a first-order phase transition. Secondly, the almost constant behavior of a{T = 0, fi) with 
increasing chemical potential is unrealistic in the absence of explicit baryon (nucleon) degrees 
of freedom including their interactions. 

What is actually required as a starting point for extensions to non-zero chemical potential 
is a realistic equation of state at finite baryon density, incorporating the known properties of 
equilibrium and compressed nuclear matter. In such a framework, the density dependence of 
the chiral condensate {ipi)) (or of the scalar field a) is well known to be quite different from the 
profile shown in Fig. 10. The magnitude of (tpijj) decreases linearly with density p [41], with a 
slope controlled by the pion-nucleon sigma term, and then stabilizes at densities above normal 
nuclear matter through a combination of two- and three-body correlations and Pauli blocking 
effects (see e.g. Ref. [42]). 

Irrespective of these comments, the nonlocal PNJL approach is obviously instructive in 



17 




Figure 6: Trace anomaly (e — 3P)/T^ shown 
as a function of temperature for the mean field 
case and with mesonic correlations added. 
For orientation, lattice data are shown from 
Ref. [12]. 



Figure 7: Comparison of the energy density 
e as a function of temperature for the mean 
field case and with mesonic correlations. 
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Figure 8: Fraction of pressure and energy den- 
sity P/e as a function of 4th root of the energy 
density e^^^. 



Figure 9: Squared sound velocity func- 
tion of 4th root of the energy density e^^^. 



modeling the chiral and deconfinement thermodynamics at = 0. Dealing with finite baryon 
density requires ultimately yet another synthesis, namely that of PNJL and in-medium chiral 
effective field theory with baryons. 



4 Summary and conclusions 

A nonlocal generalization of the two-flavor PNJL model, a synthesis of Polyakov loop dynamics 
with the Nf = 2 Nambu and Jona-Lasinio model, has been derived with the aim of identifying 
dominant quasiparticle deegrees of freedom and developing insight into the symmetry breaking 
scenario of QCD thermodynamics. This nonlocal approach has the principal advantage that it 
does not require the momentum space cutoff typical of NJL type models with local four-fcrmion 
couplings. The nonlocality of the interaction generates a momentum dependent dynamical 
quark mass, M{p). Important QCD constraints can now be directly implemented, e.g. through 
results from Landau gauge Dyson-Schwinger or lattice QCD calculations, or in contact with the 
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0.3 



Figure 10: (Normalized) chiral order parameter aja^ shown in the (T,fi) plane. For small 
chemical potentials a crossover transition is manifest. For large chemical potentials a first order 
phase transition is apparent which terminates in the critical point (CEP) at (JcePjA'cep) = 
(167MeV,175MeV). 



instanton model, and in combination with the known high- momentum limit of M{p). Remaining 
parameters are determined by reproducing pion properties in vacuum. Low-energy theorems 
based on chiral symmetry are shown to be rigorously fulfilled. The Z(3) symmetric effective 
potential governing Polyakov loop dynamics is fixed by comparison with the pure-gauge QCD 
equation of state from lattice QCD. Overall color neutrality is guaranteed by construction. 

The calculated chiral condensate and Polyakov loop as function of temperature show the 
dynamical entanglement of chiral and deconfinement transition which is characteristic of the 
PNJL approach: in the absence of quark-Polyakov loop coupling, the chiral phase transition 
(second order in the chiral limit) and the deconfinement phase transition (first order in pure 
gauge QCD) are separated in their critical temperatures by more than 150 MeV. Once the 
Polyakov loop couples to the quarks and suppresses their quasiparticle Green functions in the 
confinement region, the two transitions (with explicit Z(3) symmetry breaking by the presence 
of quarks and explicit chiral symmetry breaking by non-zero quark masses) become crossovers 
at a common transition temperature Tc — 0.2 GeV. The present nonlocal PNJL approach 
demonstrates this as a general feature, without the need of introducing an artificial momentum 
cutoff. 



Further results, such as the "interaction measure" e — 3P and the sound velocity Vs = 




have been computed both at mean field level and with inclusion of mesonic (quark-antiquark) 
modes in the pion and "sigma" channels. As expected, a substantial contribution to the pressure 
and energy density below Tc comes from propagating thermal pions. Around and above Tc the 
mesonic modes quickly disappear and the thermodynamics is dominated by quark quasiparticles, 
with momentum dependent dynamical masses, and Polyakov loop degrees of freedom. As a 
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schematic extension to non-zero r^iiark chemical potential fx, the chiral phase diagram in the 
(T, n) plane has been calculated and discussed. 

The nonlocal PNJL model thus proves to be a good starting point for modeling QCD ther- 
modynamics. It has the distinct advantage over previous (local) versions that the input (at the 
level of the dynamical quark mass function) can now be directly constrained by QCD. 

Forthcoming work will be directed towards extensions with A^f = 2-|-l quark flavors such that 
direct comparisons can be performed with corresponding realistic lattice QCD computations. 
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A Bosonization of the nonlocal Lagrangian 



This appendix briefly outlines how the bosonized form of the Euclidean action S^°^, Eq. (2.6), 
is derived from the action Se, Eq. (2.4). For this purpose we multiply 



/ 



by the constant M ^ := f ^a^n exp ( - / d^zC(z) [ d 



4^ 'ficix)<p*Jx) 



X 



2G 



(A.l) 



with the mcsonic fields 



faix) = ((t(x), vf (x)). The "shift" ipa{x) ^a{x) + G^(x + I) r'Q,V'(a; — I) allows then to 
"complete the square" in the exponential of the generating functional and to eliminate the 
quadratic term in the nonlocal currents. Note that we have to demand 7o-fa7o = -^aj as satisfied 
by the operators F considered in this work. We end up with 



/ 



d xd yV'(^) 



5^ {x - y) {-\$ + mg) + C{x - y)r^ Re ipo 



x + y 



V'(y) 



The generating functional (A.l) can now be written as 



exp s — / d^x d^y ip{x) 



{x - y) + mq)+C{x- y)ra Re 



X + y 



(A.l') 

and the integration over the fermionic degrees of freedom, ip and ip, can be performed, leading 
to 

x + y' 



Z =J\fdet 



(5^^) (x - y) {-10 + ma) + Re ipa 



C(x - y)ra 



X exp 



1 

2G 



d'zC{z) / d^x(^«(xX(x) , 



(A.l") 



where det denotes the functional determinant. 

For the calculations done in this work it is more convenient to write integrals in momentum 
space than in coordinate space. Since a Fourier transform is a Hilbcrt space isomorphism, the 
integrals may be expressed in momentum space without further changes. A straightforward 
calculation leads to our final result 



Z =M ^aWdet 



(271)^5(4) {p - p') {-f + m,) + \{<i>^(p - p') + C(P - P'))ra C f ^ 



X exp < — — 



1 f d^p 



2G J (27r)4 

From this generating functional the bosonized action can be read off: 



5^°^ = - In det 



+ 



(27r)^5W {p - p') {-p + m,) + \{ct>^{p - p') + <t>l{p- p'))r^ C f ^ 



d^ 



2G J (27r)4 
as indicated in Eq. (2.6). 



4>aip)(p*aip) 



(A.2) 



+ 



(A.3) 
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B Taylor expansion of the Euclidean action S^°^ 



In this appendix we derive the Taylor expansion for the bosonized Euchdean action S^°^ up 
to terms of second order. This corresponds to a systematic perturbative expansion and to the 
derivation of Feynman rules for the mesonic degrees of freedom according to the method of the 

"effective action". 

We start with a reminder of the Taylor expansion of a functional T = T[f{x)] about /(O) 
with a (small) fluctuation h{x). One has 

5T[/(x)] 



T[f{0) + h{x)]=T[f{0)]+ I dyi 



Sfiyi) 



/=/(o) 



h{yi)+ 



+ 



dyi dj/2 



S'T[f{x)] 



Sfiy2)Sf{yi) 



(B.l) 



/=/(o) 



h{yi)h{yi) + . . . 



This formalism is now applied to S^°^. In this case, /(O) = {a, 6) and h{x) = {Sa{x),6TT{x)). 
The following relation is useful: 

In det = Tr In , 

which holds both in functional and matrix space. We treat the functional space first, hence we 
decompose det = det Det or Tr = Tr (g) tr, where operators with a tilde act exclusively on 
functional and Det,tr solely on Dirac-, flavor- and color space. 

Then, the zeroth order term is easy to calculate since a is a Lorentz invariant and, hence. 



proportional to unity in functional space, i.e. cr{p'\p) = c-^^i^ 

15(4) ( 



[P-P 



Therefore, the func- 



tional trace simply gives an integration over p and the argument of the logarithm is given by 
Det [-:^ + mg + C{p)a] = {p^ + M'^ {p))^^f^'= . Prom this the mean field result (2.9) is deduced. 
The linear term vanishes by definition of the mean fields. 



In order to determine the term of second order we calculate 



rTr In A 



a{x)=a 



As 



5a{k) Sa{e) 

a continuation of matrix multiplication the logarithm of an operator A is treated as a power 
scries where the multiplication is given by the convolution. Therefore all matrix identities can 
be adopted in functional space, too, if all matrix products are replaced by convolutions. In 
particular, in this article we use operators O = 0[f] that are supposed to fulfill necessary 
convergence criteria so that the following identities hold: 



Sf{y) 
6 

Sf{y) 



Tr In 



Tr 



Oifix)] 



Tr 



Sf{y) 



-Tr 



o-'[ff4^d-'[f] 



Sfiy) 



(B.2) 



In the case of the operator A, Eq. (2.7), we may write 
5^ 



-Tr InA = -Tr 



(C (• - i))A-\C (• - k))A-^ 



6a{k) 6a{e) 

where the • stands for the arguments of the delta-function. Taking into account 



{p\A\p') 



(27r)^5(^) (p - p') [-f + mq + C{p)a\ , 



{a,7f)={<T,0) 

one carries out the convolutions and arrives at the result: 

X [-f + mq+ C{p)a\ . 

Next, calculate the traces over Dirac-, flavor- and color space; moreover the expansion of the 
quadratic term in cj)a is trivially done, leading altogether to the result stated in Eq. (2.10). 
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C Pion decay constant 

Here we present the full expression for the pion decay constant defined as 

(0|<i(0)|7r,(p))=i/,(5,,p^, (2.21) 

where ^ denotes the axial- vector current. 

In order to calculate the axial current matrix element one has to gauge the nonlocal action 
in Eq. (2.4). This requires not only the replacement of the partial derivative by a covariant 
derivative, 

where Tj are the Pauli matrices and A^^ {i G {!,... ,3}) are axial gauge fields, but also the 
connection of nonlocal terms through a corresponding parallel transport with a Wilson line. 



W(x,y) =:Pexp|^^ da-izTiAl{x+{y-x)a){y^j,-x^)^ , 

where we have chosen a straight line that connects the points x and y. Exemplarily, this means 
that expressions of the form 'il>{x — z) C{z) 'ip{z — y) in the action 5e, Eq. (2.4), have to be replaced 
by ^{x — z) W{x — z, z)C{z) yV{z^ z — y) il){z — y), which guarantees the (local) gauge invariance of 
the underlying Lagrangian. Following the steps described in Sect. 2, one finds that the only term 
affected by the gauging is the fermion determinant A, Eq. (2.7), which then reads in coordinate 
space: 

A^{x,y)= {^''^$y + ]p5TiJi 5{x-y)+ 

+ c(.-,)w(x,^i±l) r„^„(i!i±l) w(^i±l,; '^'^ 

where stands cither for Fq := 1 or = i^^Ti, and (fa accordingly for either the scalar field, 
£7, or a pseudoscalar field, VTj. 

The desired matrix element then follows from the gauged fermion determinant according to 



(0|JL(0)|7r,(p)) 



5^ In det A 



G 



5T,,{p)5A^^{t) 



(C.2) 

t=0 



where we have used the unrenormalized pion field, ttj = g-jiqqTTj. The complete calculation is 
rather lengthy, but straightforward and uses the operator identities (B.2) and 



5W{x,y) 



2 70 



da 75ri 5{x + (y - x)a) (y^ - x^) 

A=0 ' ■ ■ 
t=0 



5AUt) 

After a Fourier transform of the resulting expression one obtains 
(0|J^,(0)K(p)) =2itr{r„r,}Tr| f ^^^^^f^] + 

+ 2i trlrj , r,-} TV J C (q) ^i^^il ) ^ 1 ^ 

+ 2i a tr {ri,ry}x 

.Tr i[\aq,^c(q-^ ' q-.+Miqt)Miq-, 



^'^ dg2 V 2 ; {qf+M^{qt)){q-\M^{q-)) J ' 



(C.3) 
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with 



(C.4) 





Now, the decay constant can be derived from the expression (C.3) and the definition, 
Eq. (2.21), by contraction with p'*, hence 



evaluated at the pion pole = —ra^. 
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